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INVESTIGATION OF WING FLUTTER AT TRANSONIC SPEEDS FOR 
SIX SYSTEMATICALLY VARIED WING PLAN FORMS 


By George W. Jones, Jr., &nd Hugh C. DuBose 
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SUMMARY 


An investigation of the effects of systematic variations in wing 
plan form on the flutter speed at Mach numbers between 0.73 and 1.43 has 
been conducted in the 26-inch Iangley transonic blowdown tunnel. The 
angle of sweepback was varied from 0% to 60% on wings of aspect ratio 4, 
and the aspect ratio was varied from 2 to 6 on wings with 45% of sweep- 
back. The results are presented as ratios between the experimental 
flutter speed and the reference flutter speed calculated on the basis of 
incompressible two-dimensional flow. This ratio, designated as the 
flutter-speed r&tio, is plotted as & function of Mach number for the vari- 
ous wings. It is found that the flutter-speed ratio increased rapidly 
past sonic speed for sweep angles of 45° and less, indicating a favorable 
effect of Mach number. For sweepback of 609, the flutter-speed ratio was 
nearly constant throughout the Mach number range of the tests. Reducing 
the aspect ratio had a favorable effect on the flutter-speed ratio which 
was of the order of 100 percent higher for the aspect-ratio-2 wing than 
for the aspect-ratio-6 wing. This percentage difference was nearly con- 
stant throughout tbe Mach number range, indicating that the effect of Mach 
number was about the same for all aspect ratios tested. 


INTRODUCTION 


There is an urgent need for experimental ddt& dealing with the prob- 
lem of wing flutter in the transonic speed range. In order to provide a 
portion of the needed information, several flutter investigations have been 
undertaken in the langley transonic blowdown tunnel. The results of the 
first of these investigations are presented in reference l and show that 
reliable flutter data can be obtained from a slotted-throat transonic tun- 
nel. In the present investigation, the flutter characteristics of a series 
of six systematically varied high-speed wing plan forms were studied at 
transonic speeds. The purpose of the investigation was to determine the 
effects of sweepback and aspect ratio on the flutter speed for Mach num- 
bers in the vicinity of 1.0. M C" 
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The flutter tests were made at 0% angle of attack over a range of 


Mach numbers from 0.73 to 1.43. ‘The systematic plan-form variation was 
accomplished by varying the sweepback from 0° to 60° on wings with an 
aspect ratio of } and varying the aspect ratio from 2 to 6 on wings 
with a sweepback of 45°. The results of the investigation are presented 
and analyzed herein. 


SYMBOLS 


aspect ratio including body intercept 


distance in wing semichords from midchord to elastic-axis 
position measured positive rearward, ex, - l 


iri — e " Exposed half-span 
geometric aspect ratio of one wing panel, Minn giresnv ine hora 
half-chord perpendicular to quarter-chord line, ft 
half-chord perpendicular to quarter-chord line at intersection 

of quarter-chord line and wing root (except for 245 wing, 

see "Methods of Analysis"), ft 


half-chord measured streamwise at intersection of wing root and 
fuselage, ft 


wing chord measured perpendicular to quarter-chord line, ft 
first bending natural frequency, cps 
second bending natural frequency, cps 
first torsion natural frequency, cps 


uncoupled first torsion frequency relative to elastic axis, 
2 

r 
eh- (Xa/ To, ) 


structural damping coefficient 


a CPS 


structural damping coefficient for bending vibration 


structural damping coefficient for torsional vibration 
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Ve/VR 
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polar moment of inertia of wing section about elastic axis, 
slug-ft^/ft 


reduced-frequency parameter,  ab/V 

length of wing panels, outside fuselage, measured along quarter- 
chord line (except for 245 wing, see "Methods of Analysis"), 
ft 

Mach number 


mass of wing per unit length along quarter-chord line, 
slugs/ft 


dynamic pressure, lb/sq ft 
nondimensional radius of gyration of wing section about elastic 
| 2\1/2 
axis, (1,/mb 
stream velocity at flutter, fps 


stream velocity at flutter estimated for us 50, fps 


component of stream velocity at flutter, normal to quarter- 
chord line, fps 


flutter-speed r&tio 


distance of elastic axis of wing section behind leading edge, 
fraction of chord 


distance in semichords from ving-elastic-axis position to wing 
center-of-gravity position 


nondimensional coordinate along quarter-chord line (except for 
245 wing, see "Methods of Analysis"), fraction of length 1 


ratio of mass of wing to mass of & cylinder of air of diameter 
equal to chord of wing, both taken for equal length along 
quarter-chord line, m/xpb* 

taper ratio, Tip chord /Chord in plane of symmetry 

angle of sweepback of quarter-chord line, deg 


air density, slugs/cu ft 
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œ angular frequency of vibration, radians/sec 

ah angular bending frequency, radians/sec 

dày angular uncoupled torsion frequency about elastic axis, 
radians/sec 

Subscripts: 

e experimental values at start of flutter 

R calculated values based on two-dimensional incompressible-flow 


theory with account taken of mode shape and sweep (corresponds 
to subscript A in ref. 2) 


stnd. based on sea-level conditions 
APPARATUS AND TESTS 


Wind tunnel.- The 26-inch Iangley transonic blowdown tunnel is 
equipped with & slotted test section, octagonal in shape, which allows 
the tunnel to operate through the transonic speed range from subsonic 
Mach numbers up to a Mach number of approximately 1.45. A plan view of 
the tunnel, with a model installed, and a cross-sectional view of the 
octagonal test section are shown in figure 1. 


A variable and continuous regulation of the air flow is allowed by 
a set of plug valves, located between a high-pressure reservoir and the 
tunnel, and operated by a single control. A quick-operating mechanism 
closes the valves in approximately 1/2 second. 


Orifice plates of different sizes may be installed downstream of the 
test section. The orifice, when choked, permits & prescribed test-section 
Mach number to be maintained while stagnation pressure (and thus density) 
is varied from that for orifice choke up to about 75 lb/sq in. Since the 
occurrence of flutter depends on air density as well as velocity and Mach 
number, this technique permits flutter to be obtained at several Mach 
numbers on the same model by the simple process of varying the tunnel 
pressure. The tunnel air temperature varies with initial reservoir con- 
ditions and expansion in the reservoir during each run. The test-section 
velocity is therefore not uniquely defined by the Mach number. 


Support system.- The wings were mounted at OH angle of attack on a 
2-inch-diameter cylindrical sting fuselage. A fixed wing root condition 
was obtained by mounting the wing with close-fitting filler blocks and 
four 3/8-inch bolts. Figure 2 shows a flutter model mounted on the sting 
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fuselage. The fuselage nose was extended into the subsonic flow region 
of the tunnel entrance cone in order to prevent the formation of a bow 
shock wave which might reflect on the model. The support system was con- 
Sidered to form & rigid mount for the wings since the mass of the support 
System was very large when compared to the mass of & wing. 


Instrumentation.- Each model vas instrumented with wire strain gages 
located on the wing near the root and so oriented that the output of the 
gages indicated the wing bending and torsion deflections. The primary 
use of the strain gages was to provide an indication of the start of 
flutter and to obtain a record of the frequency of wing bending and 
torsion oscillations. Some estimates of the magnitude of the deflections 
and the phase angle between bending and torsion could also be obtained 
from the output of the strain gages. 


During the tests a multichannel recording oscillograph was used to 
make recordings of the strain-gage signals, tunnel stagnation pressure 
and temperature, and test-section static pressure. A sample test record 
is given in figure 3 in which the start of flutter is shown by the change 
in wing oscillations from a random form to a sine wave, the amplitude of 
which rapidly increases. 


A high-speed 16-mm motion-picture camera (approximately 1,000 frames 
per second) was used to obtain a visual record of wing deflection during 
flutter. These films served as a supplement to visual observation of the 
mode shape and magnitude of flutter. 


Tests.- The objectives of the wind-tunnel test program were to deter- 
mine the flutter speed and flutter frequency of each wing at OO angle of 
attack for several Mach numbers in the transonic range. The procedure 
followed in obtaining model flutter at a particular Mach number was to 
increase the stagnation pressure gradually until flutter was seen by an 
observer looking through a porthole in the side of the tumel. The 
stagnation pressure and Mach number were then held constant for a few 
seconds at initial flutter conditions, after which the air flow was 
quickly stopped in an effort to save the model from destruction by flutter. 


MODELS 


Model description.- Flutter tests were conducted on six high-speed 
wing plan forms which had systematic variations of sweepback and aspect 
ratio. The other geometric parameters were held constant. Four of the 
wings had an agpect ratio of 4 and angles of sweepback of the quarter- 
chord line of 0”, 45%, 52.59, and 60%. The other two wings had sweep- 
back of the quarter-chord line of 45% and aspect ratios of 2 and 6. All 
the wings had NACA 65A004 streamwise airfoil sections, a taper ratio of 
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0.6, and a ratio of body-cross-section area to wing-plan-form area of 
0.15. The model spans were 0.808 foot, 1.14 feet, and 1.40 feet for the 
wings of aspect ratio 2, 4, and 6, respectively. Drawings of the various 
plan forms tested are presented in figure 4+. Each plan form is designated 
by a three-digit number; the first digit gives the aspect ratio and the 
last two digits give the angle of sweepback. 


It was necessary to employ various materials and types of construc- 
tion in the models so that flutter could be obtained within the range of 
air densities available in the transonic blowdown tunnel. In general, 
wings having the higher length-chord ratios required stiffer structures. 
The 400 and 445 wings were made of solid compreg (a laminated, compressed, 
resin-impregnated maple). The 452 wing was made of compreg with a 
0.006-inch sheet of Fiberglas wrapped around the outside and bonded to 
the wood; the 460 wing was constructed similar to the 452 wing except 
that the thickness of the Fiberglas was 0.018 inch. In the construction 
of the 452 and 460 wings, an attempt was made to undercut the wood before 
wrapping the wings with Fiberglas; however, the airfoil shapes that 
resulted had an average maximum thickness of 5 percent of the streamwise 
chord rather than the intended 4 percent. The 245 wing had a tapered spar 
of pine 2 percent thick, with grain direction parallel to the quarter- 
chord line. This spar was sandwiched between two layers of balsa 1 percent 
thick with grain direction parallel to the airstream. The 645 wing was 
made of solid magnesium. 


Physical parameters.- Measurements were made of the following physi- 
cal parameters: elastic-axis location, first and second bending and first 
torsion natural frequencies, mass variations along the span, moment-of- 
inertia variation along the span, center-of-gravity location, and the 
structural damping coefficient in bending. A brief discussion of the 
methods used to obtain the measurements follows. 


For the determination of the elastic axis, the wings were clamped 
along a line perpendicular to the quarter-chord line and passing through 
the intersection of the wing trailing edge and the root. At several span- 
wise stations, the wings were loaded at a number of points along a line 
perpendicular to the quarter-chord line. The chordwise position of load 
application, for which no rotation of the line perpendicular to the 
quarter-chord line occurred, was determined. A straight line faired 
through these flexural centers was taken to be the elastic axis. An 
exception was the 245 wing which was clamped at the root in a streamwise 
direction. It was impractical to clamp the 245 wing in the same manner 
as the other wings because of the small length-chord ratio. This change 
in clamping caused the elastic axis to shift rearward with the result 
that the values of the elastic-axis location and the radius of gyration 
for the 245 wing were greatly different from those of the other plan 
forms. 
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The frequencies of the pertinent natural modes of each wing were 
determined with the wings clamped at the root in & streamwise direction 
which corresponded to the wing mounting during the tests. The several 
methods used to determine natural frequencies were: 


(1) Analysis of the sound response obtained by striking the wing. 
In this method, the natural frequency was obtained by striking the wing 
and tuning the sound analyzer for maximum response. The node line for a 
particular mode could be detected by the occurrence of & minimum response 
in that mode when the point of striking coincided with a point on the 
node line. 


(2) Detection of resonant frequencies with & vibrator applied at the 
wing root. The technique employed in this method involved the use of 
sand to reveal node lines and indicate the corresponding resonant fre- 
quencies of the vibrating wing. 


(5) The use of the wing strain gages to measure the frequencies of 
the decaying oscillations following release of the wing tip from a 
deflected position. 


(4) For the determination of the lower first bending frequencies 
only, the observation of vibrations under stroboscopic lighting. 


Difficulty was experienced in defining the frequency of the first 
torsional modes of both the 445 wing and the 452 wing. For the case of 
the 415 wing, three resonant torsional frequencies were found in the 
same range, two of which had equal response but node lines near the 
leading and trailing edge, respectively. The third mode had a slightly 
weaker response but a node line nearly parallel to the quarter-chord 
line and just forward-of the midchord line. For the 452 wing, a response 
of &lmost uniform amplitude over & range of frequencies was obtained. 

For both wings, the frequency corresponding to the node line most nearly 
parallel to the quarter-chord line was taken to be the first coupled 
torsion frequency. In all cases, the uncoupled torsion frequency was 
obtained from the relation given in reference 2 and included in the list 
of symbols herein. The uncoupled bending frequency was taken to be the 
same as the measured first bending frequency. It was found that, for the 
wings tested, differences in frequencies between left and right wing 
panels and between two models of the same design were smaller than the 
accuracy with which the frequency was measured. Therefore, only one 
value of the bending and torsion frequencies is presented for each wing. 


A model of each plan form was cut into strips perpendicular to the 
quarter chord and approximately 0.5 inch in width. The number of strips 
varied from 3 on the 245 wing to 11 on the 645 wing. The variation of 
mass per unit length along the span was found by weighing the strips. 
The variation of moment of inertia along the span was found by swinging 
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the strips in a calibrated torsion pendulum &bout an axis coinciding with 
the wing elastic axis. The section center-of-gravity location was deter- 
mined by balancing the strips on a knife edge held parallel to the quarter- 
chord line. 


The structural damping coefficients in bending were found by taking 
the logarithmic decrement of the bending strain-gage traces after the wing 
was depressed at the tip and released. These measurements were made in 
air. í 


Values of the parameters describing the geometric and physical prop- 
erties of the models can be found in tables I and II. The ranges of varia- 
tion of some of the more important parameters are: center-of-gravity 
oke to O.h6c, ratio of first bending to first torsion frequencies 
squared 0.0083 to 0.420, and structural damping coefficient in first 
bending 0.013 to 0.030. 


METHODS OF ANALYSIS 


The purpose of this investigation is to show the effect of changes 
in wing plan form on the flutter speed in the transonic range. Because 
of the manner in which the experimental investigation was made, however, 
the value of the mass-density parameter pu varied for the different Mach 
numbers at which flutter was obtained on the various wings. Furthermore, 
the value of the torsional frequency dh, as well as the nondimensional 
quantities Xy, 8, rg, and _ ah/a,, varied for the different wing plan 
forms. Consequently, a true indication of the effects of plan form and 
Mach number cannot be obtained merely from & comparison of the experi- 
mentally determined flutter speeds. In an effort to separate the effects 
of plan-form variations from the effects of the other variables, the 
results will be presented in the form of & ratio of the experimental 
flutter speed to a calculated flutter speed Ve/VR where VR is the 
calculated reference flutter speed and Ve is the experimental flutter 
Speed. The method for obtaining the reference flutter speed will now be 
discussed briefly. 


Reference flutter speed.- The reference flutter speed as computed in 
the analysis is based on two-dimensional incompressible-flow aerodynamic 
coefficients. The method used for calculation of VR was an application 
of that given in reference 2 (VR corresponds to VA in ref. 2). The 


terms in the analysis which involved the spanwise derivative of the 
velocity potential (bracketed terms in ref. 2) were omitted. Except in 
the case of the 245 wing, the effective wing root and tip are defined in 
the present analysis as the perpendiculars to the quarter-chord line at 
the intersections of the quarter-chord line with the actual root and tip, 
respectively. In the case of the 245 wing, the effective root is defined 
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as the perpendicular to the elastic axis at the intersection of the elastic 
exis and the root, and the effective tip is defined as the perpendicular 
to the elastic axis drawn through the intersection of the half-chord line 
and the wing tip. In all cases, the root semichord br is one-half the 
effective root chord; and the length 2 is the length of the effective 
wing panel, that is, the perpendicular distance between the effective 
root and the effective tip. 


Two modes were used in the analysis. The frequencies used were the 
measured. first bending frequency and tbe uncoupled first torsion fre- 
quency. The mode shape of the wings during flutter was represented in 
the analysis by the first bending and first torsion mode shapes of a 
uniform cantilever beam. The mode shapes of a uniform beam were thought 
to approximate the flutter mode shape with a sufficient degree of accuracy 
and were employed in preference to the mode shapes of a tapered beam as a 
matter of convenience. There were two indications, however, that higher 
netural modes than the first should perhaps have been included in the c&l- 
culations for the wings having the higher length-chord ratios. High-speed 
camera results indicated, as will be discussed later, an outward displace- 
ment toward the tip of the region of maximum curvature of the flutter mode 
Shape for these wings. This suggests the presence of higher mode shapes 
than the first, particularly in bending, in the flutter mode. Also for 
the higher length-chord-ratio wings, namely the 452, 460, and 645 wings, 
the flutter frequency was near the second natural bending frequency which 
in turn was less than the first torsion frequency. 


The values of k were weighted along the span in accordance with 

the chord variation. The spanwise variation of the Theodorsen functions 
F(k) aná G(k) were approximated by a linear variation between their 
root and tip values. With the Theodorsen functions represented in this 
manner, it was possible to set up the analysis so that computation of the 
flutter determinant coefficients and solution of the determinant could be 
accomplished quite rapidly by automatic, punchcard computing equipment. 
À comparison of the results obtained on the automatic equipment using the 
linear variation of F(k) and G(k) with those by manual methods using 
the actual values of F(k) and G(k) at 10 points along the span showed 
excellent agreement. 


The solution of the flutter stability determinant was carried out in 
such a way that a curve of structural damping coefficient g against 
Vn/Or% was obtained for each air density. Two such curves are shown in 
figure 5. The smooth curve designated "Typical 400 wing curve" is typical 
of most flutter experience. On this curve, an increase in damping from a 
value of zero results in an increase in VR since VR varies directly 
with Vn. A common practice when the curves of g are of this type is 
to present VR values based on a structural-damping value of zero which 
gives a value of VR that is conservative when compared with values 
based on actual damping. In the present analysis, a type of curve 
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frequently encountered is shown by the "S" shaped curve in figure 5 
designated "Typical 445 wing curve." Examination of this curve shows 
that an increase in damping from & zero value results in & decrease in 
Vg. If zero damping is assumed in this case, the value of Vp obtained 
is unconservative when compared with a value based on the actual damping. 
For this reason, the values of VR employed herein are based on the 
measured value of the structural damping. Since only the structural 
damping in bending was measured, it was assumed that gh = & = Es 


Flutter-speed coefficient.- In order to provide some indication of 
the manner in which the actual flutter speed might be expected to vary 
with wing plan form and Mach number for wings having & constant value 
of the density-ratio parameter pu, a flutter-speed coefficient V' / bg, 
was determined from the theoretical and experimental data. 


The parameter V'/bsa, was determined as follows: The theoretical 
value of the flutter-speed coefficient VR/ bg, corresponding to a value 
of the relative-density parameter u of 50, was determined for each of 
the wings tested from the data given in tables I and III. In order to 
obtain values of  Vg/bsaày corresponding toa y of 50 for the 245, 

460, and 400 wings, it was necessary to extrapolate the straight lines 


given by plots of VR/ ben, against Vụ. On the basis of the assumption 
that the curves of Va/Va are independent of qu, the theoretical value 
of Vn/bsœ, for each wing was corrected by use of the values of Va/ VA 


as multipliers. The products obtained by this procedure are the values 
of V'/bs%k which will be presented as a function of Mach number. It 
should be noted that bg is the half-chord of the wing taken at the 
intersection of the wing and fuselage. 


RESULTS AND DISCUSSION 


General comments.- Visual observations, examination of high-speed- 
camera results, and comparison of flutter frequencies with natural fre- 
quencies indicated that the flutter obtained in the current tests was 
most probably of the classical bending-torsion type. A sequence of 
views of flutter on two swept wings is presented in figure 6. One 
sequence shows the Jo wing viewed from a position downstream of and 
slightly below the wing; the other sequence shows the 460 wing viewed 
from a position upstream of and slightly above the wing. These edge-on 
views give some indication of the mode shape during flutter. Observation 
of the 460 wing, which has a higher length-chord ratio than the l5 wing, 
shows that the mode shape in bending is characterized by an outward dis- 
placement from the wing root of the region of maximum curvature. This 
shape suggests the presence of components of natural bending modes higher 
than the first in the flutter mode and therefore suggests that perhaps 
higher modes should be considered in the analytical representation of 


the flutter mode. 
ooo o 
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Observation and time records of strain-gage traces showed that in 
the majority of cases both wing panels of each wing fluttered almost 
simultaneously. In the few cases where simultaneous flutter did not 
occur on both wing panels, & separate point is presented in the results 
for tbe beginning of flutter on each panel. 


The results of the investigation are presented in table III and 
figures Y to 10. Table III contains the values at flutter of the 
pertinent physical quantities such as experimental flutter speed Va, 
and parameters such as mass-density ratio u and flutter-speed coeffi- 
cient Ve/b,Q,. Also given are the values of the calculated reference 
flutter speed VR. In tbe figures, VR is used as a normalizing factor 
in the graphical presentation of the experimental flutter speeds. 


Effects of sweep on flutter-speed ratio.- The effects of sweepback 
on the variation of flutter-speed ratio with Mach number for wings with 


an aspect ratio of 4 and sweepback of 0%, 45%, 52.59, and 60° are shown 
in figures 7 and 8. The curves of Ve/VR against Mach number, presented 
in figure 7, are well-defined by many data points for the O9, 52.59, and 
609 sweep angles. The curve for 45° sweep lacks points between Mach num- 
bers of 1.0 and 1.4 and, as a consequence, the fairing of the curve in 
this region is somewhat arbitrary. More data are needed for the 600 
sweep curve below Mach number 1.0. The curve is extrapolated, however, 
along a reasonable path to a Mach number of 0.8 as shown by the dotted 
line, 


There is an indication of a variation with sweep angle of the values 
of Ve/VR at Mach numbers less than 1.0. The calculated flutter speeds 
at these Mach numbers are fairly close to the experimental values for 
the wings of O9, 45°, and 52.59 sweepback; whereas a somewhat greater 
difference between theoretical and experimental values is indicated for 
the wing of 60% sweepback. A similar trend in the magnitude of the 
variation of Ve/VR with sweep at moderate subsonic Mach numbers is 
shown in reference 2, figure 15. The actual values of Ve/VR in fig- 
ure 7 are, however, lower than the values of reference 2. 


In order to show more clearly the effect of sweepback on the varia- 
tion of flutter-speed ratio with Mach number, the curves of figure 7 have 
been normalized by dividing each value of Ve/VR by the value of Ve/VR 
at a Mach number of 0.8. These normalized curves are presented in fig- 
ure 8. The curves show that the favorable increase in flutter-speed 
ratio with Mach number in the transonic and low supersonic range is less 
pronounced as the sweep angle is increased; that is, the compressibility 
effect, although first appearing at about the same Mach number for each 
wing, is consistently less for the more highly swept wings. 


. Effects of aspect ratio on flutter-speed ratio.- The effect of aspect 
ratio on the flutter-speed ratio for the wings of 45% sweepback is shown 
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in figures 9 and 10. The flutter-speed ratio Ve/Vr is plotted against 
Mach number in figure 9 for wings having aspect ratios of 2, 4, and 6. 
There is a significant decrease in Ve/VR with increasing aspect ratio 
at Mach numbers of the order of 0.8. The values of Ve/VR are 1.52, 1.0, 
and 0.77 for aspect ratios of 2, 4, and 6, respectively. This variation 
of Ve/VR with aspect ratio is in agreement with trends shown in refer- 
ence 3, particularly for the lower aspect ratio wings. It should per- 
haps be pointed out that the relatively high values of Ve/VR noted for 
the aspect-ratio-2 wing may be attributable in part to differences in 

the method of obtaining the elastic-axis position of the 245 wing as 
compared to that for the other wings. As mentioned previously, the 
manner of clamping this wing in the experimental determination of the 
elastic-axis position, resulted in & comparatively far rearward elastic- 
axis position with consequent large values of (rą)? and |x„|. A 
sample VR calculation for the 245 wing employing an assumed elastic- 
axis position and effective root chord definition similar to that of the 
other wings, with corresponding values of (ra IS and Xq, gave values 

of Ve/ VR approximately 10 percent less than the values shown in 

figure 9. The unconservative values of VR as shown in figure 9 for the 
wing of aspect ratio 6 are disturbingly high and, as yet, no satisfactory 
explanation can be given for the degree of unconservatism. One possi- 
bility is that the inclusion of higher modes in the analysis than first 
bending and first torsion may lower the values of Vp. Some thought has 
also been given to the possibility that the flutter encountered on the 
645 wing was of the single-degree-of-freedom bending type described in 
references 4 and 5. This possibility seems unlikely, however, on the 
basis of the curves given in figure 3 of reference }. 


The curves of figure 10 are the curves of figure 9 normalized with 
the value of Ve/VR at a Mach number of 0.8. Since the curves nearly 
coincide, they indicate that the variation with Mach number of the flutter- 
Speed ratio in the transonic and low supersonic range is practically the 
same for aspect ratios 2, 4, and 6. 


Application of results.- In this section an application is made of 
the curves of Ve / VR to obtain curves of flutter-speed coefficient, 
V'/bs%, against Mach number for a given value of mass-density ratio y. 
The curves of V'/bga, as a function of Me, which are presented in fig- 
ures ll and 12, give an indication of the effect of wing-plan-form 
geometry on the variation of the flutter-speed coefficient with Mach 
number. Some of the differences between the curves for the various wing 
plan forms is also attributable to variations in certain of the physical 
parameters which describe the wings. Perhaps the most significant differ- 
ence in the physical properties of the wings is found in the bending- 
torsion frequency ratio which varies from about 0.36 to about 0.098 as 
the sweepback of the aspect-ratio- wings is varied from 0% to 600, and 
which varies from about 0.65 to about 0.09 as the aspect ratio of the 
150 sweptback wing varies from 2 to 6. The comparison of wing plan forms 
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shown in figures 11 and 12 is, therefore, one in which changes in fre- 
quency ratio accompany changes in wing-plan-form geometry. The direction 
in which the frequency ratio changes with plan-form geometry for the 
curves shown in figures ll and 12 is perhaps the same, however, as might 
be expected on actual airplane wings without external appendages. 


Comparison of the curves of figures ll and 12 indicates that for a 
given root chord and torsional frequency there is a relatively small 
effect of sweep on the flutter speed for Mach numbers less than 1.0. At 
higher Mach numbers, however, there is & large systematic effect of 
sweepback on the flutter speed with the straight wing having the highest 
flutter speed. The curves of figure 12 show that for the same root chord 
and torsional frequency, there is a very large increase in flutter speed 
as the aspect ratio is decreased from 4 to 2, but that a relatively small 
increase in flutter speed accompanies & reduction in aspect ratio from 6 
to 4. "These comparisons are made on the basis of a fixed torsional fre- 
quency and chord length as the plan form is varied. In practice, however, 
it would seem likely that the torsion frequency would vary vith plan form. 
In any case, the comparative values of the actual flutter speed for the 
wings of different plan form will obviously vary with the torsion fre- 
quency. Again, the differences in the method of obtaining the elastic- 
exis position for the 245 wing should perhaps be considered in comparing 
the values of flutter-speed coefficient for the 245 wing with the values 
for the other wings tested. In the case of the flutter-speed coefficient 
V'/bgt, the effect of the rearward location of the elastic axis obtained 
for the 245 wing is evident in the low values of uy resulting from the 
use of the relatively large values of Lal in the formula for the 
uncoupled frequency. 


It should be emphasized that the comparisons of the effect of plan 
form on the flutter speed as shown in figures ll and 12 are based on wings 
which are characterized by particular values of the various pertinent 
physical parameters such as frequency ratio and center-of-gravity position. 
Changes in the values of these parameters for one wing with respect to 
another would, of course, alter these comparisons. For example, a rear- 
ward movement of the center-of-gravity position of the straight wing would 
lower the flutter speed of this wing with respect to the flutter speeds 
of the swept wings. 


The manner in which the curves of figures 11 and 12 were obtained 
illustrates the application of the curves of Ve/ Vp given in figures T 
and 9. Caution should be exercised, however, in &pplying these flutter- 
Speed ratios to the determination of the flutter speed of wings which 


have values of 0/0, Xas 8, r,, and p much different from those 


which characterize the wings of the present investigation. It might be 
hoped that the reference-flutter-speed calculations, &s obtained in the 
present paper, have adequately removed from the results the effects of 
such variables as the center-of-gravity position, and that the curves of 
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V./Vp against Mach number are a function of plan form only. It is not 
entirely evident, however, that such is the case and it is thought that 
further investigation of particular wing plan forms having different 
values of the various pertinent parameters which go into the reference- 
Speed calculation &re required in order to establish the applicability 
of the results obtained. 


CONCLUSIONS 


The results of & systematic study of the effects of variation in 
wing plan form on flutter at transonic speeds indicated: 


l. The flutter-speed ratio increased rapidly past sonic speed for 
sweep angles of 45° and less, indicating a favorable effect of Mach num- 
ber. For sweepback of 60%, the flutter-speed ratio was nearly constant 
throughout the Mach number range of the tests. 


2. Reducing the aspect ratio had a favorable effect on the flutter- 
speed ratio which was of the order of 100 percent higher for the aspect- 
ratio-2 wing than for the aspect-ratio-6 wing. This percentage difference 
was nearly constant throughout the Mach number range, indicating that the 
effect of Mach number was about the same for all aspect ratios tested. 


5. Further investigation is needed to establish the generality of 
the results employing the reference flutter speed as a normalizing factor. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 18, 1957. 
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TABLE III.- EXPERIMENTAL AND ANALYTICAL RESULTS 
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TABLE III. - Continued 
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TABLE III.- Concluded 
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Figure 1.- Plan view of Langley transonic blowdown tunnel with flutter 
model installed. 
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Figure 2.- Flutter model of a wing with an aspect ratio of 4 and 60° of 
sweepback, mounted on sting fuselage. 
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Figure 3.- Sample oscillograph record of flutter test. (Wing 445 at „ 
M = 0.813.) 
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Figure 4.- Plan forms of flutter models giving aspect ratio, sweep angle, 


and model designations. 
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Figure 5.- Typical curves of structural damping g against the 
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445 wing 460 wing | L-80289 
Figure 6.- Sequences from high-speed motion pictures of flutter 
of 445 and 460 wings. 
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Figure 7.- Effect of sweepback on variation of flutter-speed ratio with 
Mach number for aspect-ratio-4 wings. 
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-4 wings. 
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Figure 8.- Effect of sweepback on variation of normalized flutter-speed 
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Figure 9.- Effect of aspect ratio on variati 
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Figure lO.- Effect of aspect ratio on variation of normalized flutter- 
speed ratio with Mach number for wings with 45° sweepback 
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Figure 11.- Effect of sweepback on variation of estimated flutter-speed 
coefficient with Mach number for aspect-ratio-! wings. 
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Figure 12.- Effect of aspect ratio on variatio 


n of flutter-speed 


coefficient with Mach number for wings with 45° sweepback. 


u = 50. 


NACA-Langley - 8-13-53 - 325 


